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Abstract. We present the modeling of partial frequency redistribution (PRD) effects for the fluo- 
rescent emission lines of molecular hydrogen, the general computational approximations, and the 
applications to planetary atmospheres, as well as interstellar medium. Our model is applied to FUSE 
observations of Jupiter, Saturn, and reflection nebulae, allowing an independent confirmation of the 
H2 abundance and the structure of planetary atmospheres. 
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PRD FLUORESCENCE MODEL 

Frequency Redistribution. FUSE observations of planetary atmospheres and re- 
flection nebulae have revealed the need to include PRD in H2 fluorescence radiative 
transfer models. The angle- averaged laboratory frame redistribution function describes 
the conditional probability that a photon absorbed at x,- Doppler widths from the center 
of line is emitted at Xf Doppler widths from the center of line / [1]. Complete 
redistribution (CRD) refers to the photons being re-emitted in the line core according 
to the Voigt profile, while PRD takes into account the changes in the line profile due to 
coherent scattering in the line wings. PRD becomes important for fast transitions and 
integrated optical depths larger than 10^. Observable effects include shifts in the peak 
wavelength in emission due to variations of the exciting spectrum over the absorbing 
line profile, a decrease in the line-to-continuum contrast, and a decrease in the number 
of photons scattered in subordinate lines (cross redistribution, XRD). 

Radiative Transfer The PRD radiative transfer problem is more complex than the 
CRD case, due to the heavy couplings in frequency in addition to spatial correlations 
and detailed balance [0,|3l]. Even restricting the geometry to the plane-parallel case, and 
assuming fixed level populations, the problem remains computationally expensive due 
to intrinsically large optical depths and fine frequency grids. We investigated a set of 
numerical methods (see Table [l]) that would allow the full treatment of the H2 molecule 
(10^ transitions) to become feasible. We find that the multilayer approximation is the best 
approach, with a computing time independent of optical depth, and no matrix inversions 
required. This method is a layer-by-layer extension of the optically thin single layer 
solution of Liu and Dalgarno [4]. The results are consistently within 10% of the Feautrier 
lambda iteration [1] for grids of at least 10 points per dex in optical depth. 



TABLE 1. Characteristics of the radiateve transfer methods used. 



Lambda Iteration 
(Feautrier solution) 



Multilayer 



Single 
layer 



Resources 
Consistency 
Convergence 
Spatial variations 



High 

Can be unstable 
Scales with optical depth 
Yes 



Moderate 
Constrained 
2-pass layer-by-layer 
Yes 



Low 
Constrained 
Single step 
No 



RESULTS 



Planetary Atmospheres. The PRD effects on resonance lines and overlapping transi- 
tions in planetary atmospheres have been discussed previously [l5|,|6|]. The current FU SE 
observations represent the first account of PRD effects in subordinate molecular fluores- 
cent lines. The lines in the Lyman (6 — v") progression of molecular hydrogen pumped 
by the solar Ly/3 show a broad asymmetric profile, consistent with the shape of the over- 
lap between the (6 — 0) P(l) line and Ly/3, and emphasize the importance of coherent 
scattering in the line wings. The FU SE spectra of the Jupiter limb and Saturn disk shown 
in Figure [Hand Figure [21 respectively, are well-fitted by the XRD model (red). For com- 
parison, the CRD model is also shown in blue. The synthetic spectra have been obtained 
by integrating the full atmospheric models, and using the multilayer approximation. 
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FIGURE 1. FUSE MORS Jupiter data, FWHM 0.08 A (black), PRD (red), and CRD (blue) models. 
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FIGURE 2. FUSE LWRS Saturn data, FWHM 0. 12 A (black), PRD (red), and CRD (blue) models. 



NGC 2023. In spite of predictions based on the infrared H2 observations ff], high 
excitation non-thermal H2 absorption detected in the UV [8], and UV H2 emission at 
1575 — 1608A observed by HUT, the molecular hydrogen lines in the FU SE bandpass do 
not match the CRD model predictions. Here we test the assumption that this discrepancy 
is due to the effects of PRD and XRD. We employ a toy model, with an incident 22,000 K 
blackbody spectrum normalized to 5000 times Harbig galactic mean [8], illuminating 
a uniform H2 cloud with rotational temperature of 1500 K, Doppler b parameter of 
1.8 km s^\ and an integrated column of 5x10^*^ cm^^. The radiative transfer was 
performed in a single layer approximation, to minimize computational resources. The 
results shown in Figure [3] represent the first XRD calculation of a fluorescent molecular 
spectrum. This exercise shows that the PRD effects will decrease the line contrast, 
making the H2 features less prominent at high spectral resolution. In future work, we 
expect stronger peak suppression in a multilayer treatment due to the effects of multiple 
scattering. 
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FIGURE 3. Comparison between the H2 spectrum predicted by the CRD (black) and PRD (red) models. 
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